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A Power-flow Based Standardized Modeling and 
Simulation Method for Complex System 

Li-Rong Wang* , Jia>Cai Wang'* and Ichiro Hagiwara* 
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MeguFo-ku, Tokyo 

^Beijing Institute of Technology, Mechatronics Research Center, School of Mechanical En^eering and 
Automation, Beijing^ 100081 , R R. China 

Abstract This paper presents a power-flow based standardized modeling method and automatic 
simulation procedure for complex engineering .system conqiosed of mechanical, electrical and * 
hydraulic subsystems. Three kinds of conventional modeling methods, bond-graph, 
fimctional-model and block-diagram, are briefly summarized first Then k standardized modeling 
and simulation method (SMSM) is presented and implemented. Due to take advantages of 
bond-grs^h, fimctional-model and block-diagram in Sunulink of commennal code of Matlab, 
SMSM blends of standardized element units, effort-flow based graphic fiame, routinized modeling 
procedure and automatic simulation without need to writing out system mathematic equation. 
Finally, three types of mechanical systems with typical linear and nonUnear performances are 
studied. Comparisons of modeling procedure and graphic models of SMSM, bond-graph method 
and &nctional-model method clearly demonstrate .differences and relationships among those 
methods. Simulational results of the SMSM models, which meets well with the numerical results of 
their matiiematic equations and fimctional-models, verify die practicability and effectiveness of 
SMSM. SMSM to conveniently create highly readable, imified and hierarchical models for 
multidisciplinary systems and to efficiently carry out automatic simulation process provides 
engineers with an effective modeling and performance prediction technique for large-scale product 
design and development 

Key words: Modeling, Simulation, Power-flow, Ftmctional model. Block diagram. Bond graph 

1 Introduction 

Nowadays engineering systems are becoming increasingly complex to accomplish more demanding 
tasks in more efficient ways, which usually contains mechanical, electrical and hydraulic 
subsystems that interact in complex ways* In design and development of such multidisciplinary 
system, prediction and analysis of its performance are very important to ensure design quality, 
shorten product design cycles, reduce development cost and explore optional schemes especially in 
concept design stage. With increasing complication of engineering system, CAB (Computer Aided 
Engineering) as a powerfiil tool to investigate performance of complicated systems becomes 
particularly necessary. Umfied modeling approach to conveniently obtain legible and hierarclucal 
models of multidisciplinary systems is the primary to CAE system technology, and automatic 
derivation of system mathematic equation and automatic system performance simulation largety 
improve efficiency of CAE technology. 

Unified modeling method and automatic numerical simulation have been studied for many years. 
Gmphic modeling method, as abridge between physics and mathematics descriptions of engineering 



system* is one solution to get a unified model legible for describing multidisciplinary systems in a 
standard way and to enable engineei^ in different fields to understand the standardized model 
without knowing deep theories in other fields. Power-flow based graphic modeling approaches such 
as bond-graph method developed by Ronald C. Rosenberg and Dean C. Kamopp [1-7] and 
functional-model proposed by Sumida Shi2Uo,Nagainata Akio [8-19] have the capabilities to obtain 
unified model for multidisciplinary system and to automatically generate mathematic equation 
directly from gmphic model. Block-diagram oriented tools of Simulitik in the commercial code of 
Mafiab [20], Auto-code Generator in the commercial code of Beacon and SystemBuilder in fhe 
commercial code of MATRDCx express mathematic equation m a graphic way and perform 
automatic simulation, which are especiaUy.useful in design and build of control system. 

In this paper a practical and reliable Standardized Modeling and Simulation Method (SMSM in 
short) based on power-fiow is proposed, which takes advantages of bond-graph, functional-model 
and block-diagram model in Simulink to establish standardized model and accon^lish automatic 
simulation for conqiUcated system. As applications of this modeling technique, several typical linear 
and nonlinear niechanical systems are numerically simulated. The analysis results as well as the 
effectiveness and ef[iciency of SMSM are verified. This method can be applied into some other 
systems such as electrical, hydraulic and their mixed systems, and also can be conabined wi^ other 
n^ndelmg and design technologies such as finite element method, experimental statistics method to 
erihance capability of CAE system technology. 

2 Comparisons of Conventional Modeling Methods 

Dynamics of physical system stems fi:dm power exchange among its components described as 
power-fiow. Power statement of the product of a pair of effort and flow in various systems shows 
striking analogies, as shown in Table I. Bond-graph and functional-model are developed on the basis 
of the analogy of power descriptions and power conservation, in which efifort and flow are used as 
state variables to connect various systems into a unified form. In this section three kinds of modeling 



Table 1 Standardized bond-variables and elements in bond-graph 





Mechanical system 


Electrical system 


Hydraulic system 




Translatoiy system 


Rotational system 


Effort 


Force: F(N) 


Toique: T (N-m) 
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Flow 
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mel]u)ds» bond-graph, fbnctioiuil-inodel and block-diagiam in Simulink, are reviewed and 
compared. 

(1) Bond-graph modeling method 

Bond-graph is a highly organized schematic modeling method using generalized symbols or 
elements to represent abnost all applicable types of lumped-parameter systems. Pictorial model of 
bond-graph, for example, the bond-gii^h model of Fig. 4 (b), is composed of line called bond and 
nine basic bond-graph elements grouped into four groups as shown in Table 1: the dissipation field 
of Resistor element; the storage fields of ^ertia element and Capacitor element; the general 
junction-stmctore groups of 0-junction and 1-junction, transformer {TF) and gyrator (GY); the 
source fields of active ports (Se and Sfj. 

Bond connecting element with junction marks a power connection with bond-variables of effort and 
flow written in the upper and lower of the bond and a half arrow pointing in the direction of power 
flow. A single mark of causal stroke on a bond indicates input-output relation between effort and 
flow at each port, which flow in opposite directions. Each passive 1-port element o^R, C and / 
represents a constitutive law between bond-variables of efifort and flow or their time mtegrals. R 
dissipates power from system, and C and / storage power. l*junction also named as common flow 
junction has the special nature that tiie flows on the bonds attached to it are equal and the algebiaic 
sum of all efiforts attached to h is zero. On the contrary, 0-junction namely common effort junction 
possesses the fecial nature of that the efforts on the bonds attached to it are equal and algebraic sum 
of all flows attached to it is zero. 0-junction is regiaided as the dual junction of 1-junction in the sense 
that the roles of effort and flow are reversed. Ideal 2-port element of 2F relates effort to effort and 
flow to flow by a proportional dimensionless coefficient called modulus of the tran$fonner. While 
ideal 2-port element of relates the efiSort at one port to the flow at the other by a dimensionaf 
parameter of effort divided by flow. All Junction 0,1, IF and GF are power conservation. Sf and Se 
supply power into Bond-graph. 

With these generalized elements and their built-in analogy in different systems, bond-graph 
modeling process can be performed in an orderly way .as shown in Table 2, which leads to economy 
of thought and helps engineers to rapidly establish an organized dynamic analysis model for 



Table 2 Bond-graph modeling procedure 



Mechanical system 


1 . Simplify physical system into a lumped-parameter model. 

2. Establish one 1-junction for each velocity. 

3. Connect /to l<juncions and insert 0-jimction between a neighboring 1-junction pair. 

4. Comiect C and R to 0-juncions. 

5 . Connect 0 and l-junctions directly or using TF or GF if necessaiy. 

6. Insert force and velocity sources appropriately. 

7. Assign half arrows and causal strolce to bond. 

8 . Simplify the graph. 


Electrical syslem 


1 .Establish' circuit construction. 

2. Establish a 0-jimction for each node with a distinct voltage in the circuit 

3. Insert each 1-port element between a neighboring pair of 0-junctions by attaching the 
1-port to a 1-junction and then bonding the 1-junction to two neigihboring O-junctions. 

4. Assign half arrows and causal stroke to bond. 

5. Choose a particular node to be ground and simpli& the graph. 


Hydraulic system 


1 . Simplify fhud line into a distributed-paiameter system. 

2. Estsdilish a 0-junction for each pressure. 

3 . Attach C element to 0-junction. 

4. Attach J, Rot IF \o 1-junction and connect the 1-junctions between two neighborins 

O-junctions. 

5. Simplify the graph. 



3 



complicated system. Another remarlcable feature of boad-graph is that standard system state-space 
equation can be automatically derived in an orderly way from a given bond-graph diagram for 
subsequent numerical analysis. 
(2) Functional-model method 

Functional-model method composed of element unit and operation is another kind of power-flow 
based modeling method, wluch demonstrates physical function of system and power-flow 
transmitted in system in a visual way and enables engineers in various specialties to easily 
understand and collaborate with each other [8]. State variables of flow and potential, whose product 
becomes power as that of flow and effort ux bond-graph> construct mam graphic fiame as shown in 
Fig.2 (b). \nsualized element units in various systems as shown in Table 3 explicitly explain the 
basic constitutive law between flow and potential. Expansion and unification of functional-models 
make it easy to set up hierarchical model from small-scale to large-scale subsystem, to assemble 
machine from parts to total product or disassemble machine and exchange parts [12]. Functional 
model can cooperate with mechanism model to control parameter and modeling nonlinear physical 
property [13]. Furthermore, system state-space equation can be derived automatically from diagram 
of fimctional-model, and system perfonnance simulation can be performed by numerical 
computation of the mathematic equation. 

(3) Block-diagram method 

The block-diagram model in Simulink of commercial code Matlab is reviewed in this section. 
Simulink block-diagram is a pictorial model to represent a known mathematical model of a dynamic 
system, which consists of Erected lines and blocks uxterconnected by lines. Each line represents a 
single variable flow and connects between block output and block input. Each type of block 
represents a particular relationship between its input and output variables. A block-diagram can 
contain any number of instances of any type of block. A simple way to create a block diagram is to 
start with a system mathematic equation. Nonlinear block library helps, to establish nonlinear modeL 
The capabilities to group blocks into subsystem and to create new block by using ^-Function allow a 
hierarchical model to be built using both top-down and bottom-up approaches for a complex system 
with complicated linear and nonlinear performances. Moreover, subsequent simulation on block 
diagram can be automatically performed by Simulink solver using numrarical computation methods 
in Matlab without writing out system eqiiation implied in block diagram [20]. 

(4) Comparison of modeling methods 

Unified modeling approadies of bond-graph and ftmctional-model are visualized descriptions for 
lumped-paramet^ physical model on the basis of power-flow, whidi have the similar advantages 
as follows: 

• General symbols are used to model multidisciplinary system in a unified format, and graphic 
models involved different energy type demonstrate striking analo^es. 

• Physical function and power interchange between components of a system are illustrated in a 
pictorial way. 

• Standard fomi of state-space equation can be automatically derived from graphical model' 

Due to these similarities, a bond-graph diagram can be routinely converted into an equivalent 
functional-model by replaciiig each bond-graph element into equivalent functional-model unit and 
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Table 3 Element unit in various dynamic system 
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operation according to the following transformation .rule: ) 
L Each bond in bond-graph corresponds to a pair of directed signals of flow and effort in 
fimctional-modeL 

2. 0-junction in bond-gmph denotes a pair of operations in fimctionaUmodel that the flows are 
summed and effort is distributed. And l-junction is treated in similar way that efforts are summed 
and flow is distributed. 

3. Causal mark of one bond in bond-graph indicates tiie directions of corresponding fLo^ and 
potential in functional-model. 

4. TF is replaced by proportionally A ^ Block diagram 
connecting two pairs of effort and 
flow wi& a coefBd^t respectively. 
G7 is replaced by relating the effort 
at one port to the flow at the other 
with a coefficient respectively. 

5. Effort source in bond-graph 
corresponds to an input potential 
signal in &nctional-model» and flow 
source in bond-graph conesponds to 
an input flow signal in 
fimctiional-model. 




SMSM 

Fig. 1 Relationship among various modeling and simulation methods 



Table 4 SMSM modeling procedure 



Mechanical system 


1. Simplify physical system into lumped-parameter model. 

2. Establish a pair of operations offeree sum and velodty distribution for each velocity. 
Connect element unit / to its proper force and vdod^ pair. 

3. Insert a pair of operations of velocity sum and force distribution between the two 
neighboring pairs of force sum and velocily distribution if there are C and jR. Then, 
connect C and jR to this pair of operation properly. 

4. Connect two pairs of effort and velocity together by multiplying a coefficient 
respectively if necessary. 

5. Establish nonlinear blodc using ^-Function and subsystem tools. 

6. Insert force and/or velocity source signal generators appropriately. 

7. Simpli^ the eraph by erouping some blocking into subsystem. 


Electrical system 


1. Establish circuit construction. 

2. Establish a pair of operations of current sum and voltage distribution for each node 
with a distinct voltage in the circuit 

3. Insert a pair of operations of voltage sum and current distribution between two 
neighboring pairs of currents sum and voltage distribution. Cormect element unit to 
this the pair of operation. 

4. Insert voltage and/or current source signal generators appropriately. 

5. Simplify the eraph by choosing a particular voltage node to be ground. 


Hydraulic system 


1 . Simplify fluid line into distributed-parameter systeriL 

2. Establish a pair of operations of volume flow rate smn and pressure disfaribution for 
each pressure. . 

3. Attach C to the operation of volume flow rate sum and pressure distribution. 

4. Insert a pair of operation of pressure sum and volume flow rate distribution between 
two neighboring pairs of volume flow rate sum and pressure distribution. 

5. Attach element units of / and it to the pair of operations of pressure sum and volume 
flow rate distribution. 

6. Connect two neighboring pairs of volume flow rate and pressure together by 
multiplying a coefficient respectively if necessary. 

7. Insert pressure and/or volume flow rate source signal generators appropriately. 



Even though bond-graph and &nctional-model have similarity and can be converted each other, they 
have their own properties. Bond-graph requires special learning due to its special symbols, which 
make its graphical models are not easily accepted by engineer. On ^e contrary, fimctional-model 
uses ordmary engineering symbols, so it can be more easily to be understood and be* applied in 
engineering. Some nonlinear problems are difQcult to be described by bond-graph, but they can be 
convenioitly considered In functional-model by adding mechanismrmodel. 
Comparmg with bond-graph and fimctional-model made up of element units, block-diagram model 
in Simulink built by mathematic synabols has much difficulty in modeling establishment directly 
j&om physical system. Thus, block-diagram describing mathematical model in a clear way is hard to ' 
demonstrate physical function and power*flow of system. 

3 Standardized modeling and simulation method 

SMSM is developed by taking advantages of bond-graph, functional-model and block-diagram as 
shown in Fig. 1. The standardized element units are subsystem blocks built by Simulink library 
blocks as shown in Ted>le 3, which adopt the concepts of bond-gr^h element and functional-model 
element symbols to represent element functional characteristics in various systems between state 
variables of effort and flow in a simple way. The proposed standardized modeling procedure 
directly modeling physical system is shown in Table 4, which combines the modeling procedure of 
bond-graph with the transformation rule convertuig bond-graph mto &nctional-model. The main 
graphic frame of SMSM model is constructed by flow and effort as flow and potential in 
functional-model, and can be built in hierarchical way by creating subsystem block. Nonlinear 
performance can also be conveniently modeled by creating nonlinear subsystem like the mechanism 
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model in fiinctional-method. Hence, SMSM is a power-flow based modeling approach, which helps 
to conveniently build up understandable and unified model for multidisciplinary system and 
provides insight into how a model is organized and how its parts interact Automatic numerical 
simulation can be implemented on SMSM gcapUcal model by Simulink solver with no need to write 
out system equation and make program to calculate the equation. Therefore, users are able to 
concentrate on modeling decisions without bemg saddled with algebraic drudgery and burden for 
making numerical computation program. 

SMSM has highly flexibiUty. It can be established directly from physical model, or converted from 
functional-model and bond graph model. Block-diagram set up from matfaematic equation or 
physical .models such as electrical circuit can be built into hierarchical subsystem as a part of SMSM 
model. Model change can also be conveniently performed on the existed graphical model. Moreover, 
SMSM simulation process is readily incorporated with other analysis and design metiiods of 
dynamic system by using Matiab toolboxes, for example experimental analysis method, 
optinuzation design method. 

To sOm MP, SMSM collects tlte advantages of highly readability, quick establishment, explicit 
frmction and power-flow description, wide applicability and automatic simulation for numerical 
study of complex systems. 



(a) Structural model 




(b) Functional-model (c) SMSM-modd 



4 Applications 

Several typical mechanical systems are studied to demonstrate SMSM modeUng and simulation 
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(b) Simulation results of SMSM-modd 
Fig.3 Simulational results ofSDOF mechanical system 

procedure. C6nq)arisons of functional-model, bond-graph model and SMSM-model are also 
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presented. These examples are kept simple enough to show abilities of the proposed approach in 
solving linear and nonlinear problems. Validity and practicality are illustrated by comparisons of 
simulational results. 

(1) A SDOF (Single Degree of Freedom) mechanical system 

As shown in Fig. 2 (a), this mechanical system is composed of a mass, a spring and two dashpots. 
The functional-model [8] and SMSM-model shown in Fig.2 (b) and (c) are similar. It is obvious that 
the SMSM-model can illustrate fimction of each element unit and direction of power-flow a? the 
fimctional-model does. However, the SMSM graphical description is more readible and sux^le than 
that of the fimctional-model The SMSM-model can be established by converting the 
functional-model accordmg to the transformation rule introduced in the last section or directly set up 
fiom the system mechanism model according to the standard modeling procedure shown m Table 4. 
The simulation results shown in Fig. 3 (b) under the force and velocity sources shown in Fig. 3. (a) 
meet with the numerical results of system stat&*space equation. Simulation of SMSM-model can be 
performed more conventientiy by simulink solver than making numerical coiz^utational program as 




(c) Fimctional-inodel (d) SMSM-model 

Fig.4 Models of suspraded body 



(2) Suspended body 

A typical vibration problem show in Fig. 4 (a) is a suspended body, which is assumed to vibrate in 
small angular and translational motion in the plane with neglecting horizontal motion. Its various 
models are shown in Fig. 4 (b), (c) and (d). 

SMSM-model shown in Fig. 4 can be tranformated from bond-graph model of Fig. 4 (c) and 
functional-model of Fig. 4 (b). It can also be directly established from the structural model of Fig. 3 
(a) as the following steps according to the modeling procedure shown in Table 4. Step 1: establish 

four pairs ofvelocity-forc6fi:ames,v^ and O) for the rigid body and andx^ for the suspension; 

Step 2: attach element units of liT, C and Af to their corresponding velolcily-force frams; Step 3: 
connect the four velocity-force frames together by sum and distribution operations and 
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transfonnation coeSicient blocks according the relationships among them: the transformation 
coef&cience blocks of X / and L2 as shown in Fig. 4 (d) connect translatary subsystem with rotational 
subsystem in the same way as TF shown in Fig. 4 (b), and the operations of sum flow of -v/ and -y? 
and effort distribution ofy} and fz m Fig. 4 (d) is eqivalent to the 0-juaction in bond gnqph of Fig. 4 
(b) respectively. The simulation results of SMSM model under the force source of are 
demonstrated in Fig. 5 (a), which is agreeable with the simulation results of fiinctional-model shown 
in Fig. 5 (b) in Ref. [9] and the numerical results of system state-space equation. 

(3) A nonlinear mechanical system 

A nonlinear mechanical system with typical material and friction nonUnearittes as shown in Fig. 6 
(a) are investigated. In SMSM-model shown in Fig. 6 (c), the two-stage spring atiffiiess is modeled 
by a subsystem JiT, in which a ^-function of Varstiffl is programmed to control the sti&ess property 
sho\m in the mechanism model of Fig. 6 (b). Friction between Mand ground is divided into two 
stages, one is static-fiiction when M doesn't begin to move and the other is coulomb 
(fynamic-friction proportional to the press upon the slide sur&ce when Af slips on ground. Nonlinear 
friction is modeled by another subsystem of friction control shown in Fig. 6 (c). Power-flow 
transformation through lever is described as two transformation coefficient blocks of R multiplied to 
velocity and force respectively. The simulation results under the force source of i*i are demonstrated 
in Fig. 7 (a), which is agreeable with the simuiationai results of fimctional-m^ 
in Ref [10] and the numerical results of system state-space equation. Thus, it is very convenient to 
built nonlinear model and cany out numerical simulation by using tools m Smmlink. 
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5 ConclusioB 

Unified modeling approach is important to performance simulation for complex system. In this 
paper, a comprehensive modeling and simulation method of SMSM based on power-flow is 
presented. This method is time saving and convenient for modeling a complicated systbm due to its 
advan&ges of standardized element unit, efiEbrt-flow based graphic firame, routinized modeling 
procedure and automatic simulation process. From the simulation results of several mechanical 
systems with typical linear and nonlinear properties, the effectiveness and practicability of the 
proposed method are verified. The success of this approach is that it takes advantages of bond-graph, 
fimctional-modei and block-diagram in Simulink by blend of modelmg procedure and automatic 
computation process. The proposed modeling techniques and simulation analysis method allow to 
creative a highly readable, unified and hierarchical model in a orderly way for a multidisciplinaiy 
engineering system composed of mechanical, electrical and hydraulic systems interacting in 
consplex ways and efEectively perform automatic simulation analysis on computer: The presented 



method can also conveniently combine with other types of analysis methods such as finite element 
analysis method and product design methods to enhance capability of CAE system technology for 
complicated engineering problems. 
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